Role of the cr-resonance in determining the convergence of chiral perturbation theory 
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The dimensionless parameter £ = M%/(16-k 2 F%), where F n is the pion decay constant and M n is the 
pion mass, is expected to control the convergence of chiral perturbation theory applicable to QCD. Here we 
demonstrate that a strongly coupled lattice gauge theory model with the same symmetries as two-flavor QCD 
but with a much lighter cr-resonance is different. Our model allows us to study efficiently the convergence of 
chiral perturbation theory as a function of £. We first confirm that the leading low energy constants appearing in 
the chiral Lagrangian are the same when calculated from the p-regime and the e-regime as expected. However, 
£ < 0.002 is necessary before 1-loop chiral perturbation theory predicts the data within 1%. For £ > 0.0035 
the data begin to deviate dramatically from 1-loop chiral perturbation theory predictions. We argue that this 
qualitative change is due to the presence of a light cr-resonance in our model. Our findings may be useful for 
lattice QCD studies. 

PACS numbers: 1 1.15.Ha,l 1.15.Me,12.38.Gc,12.39.Fe 



Chiral perturbation theory has been successful in explaining 
a variety of experiments involving low energy pions ( 1] . It is a 
low energy effective field theory that captures the chiral sym- 
metry properties of QCD. The dynamical properties of QCD 
are encoded through a series of low energy constants. At the 
leading order there are two low energy constants: F the pion 
decay constant and X the chiral condensate, both evaluated in 
the chiral limit. One of the important topics of research today 
is to compute these and other higher order low energy con- 
stants from first principles using lattice QCD J2, S H • In- 
terestingly, the effects of a small quark mass m, which breaks 
the chiral symmetry explicitly, can also be taken into account 
and physical quantities can be expressed as a power series in a 
dimensionless parameter £ = M 2 / '(16-7T 2 F%) where M x is the 
physical pion mass and F v is the physical pion decay constant. 
This series, which we refer to as the "chiral expansion", not 
only contains powers of £ but also powers of £ log £, £ 2 log £ 
and so on. In QCD we can estimate £ ~ 0.015 assuming 
~ 140 MeV and F^ ~ 90MeV. For later convenience we 
also define f = M 2 /(16tt 2 F 2 ) where M 2 = mE/F 2 is the 
square of the pion mass to the leading order in the quark mass. 
It is easily verified that £ w £' + 0(£' 2 ), which means that to 
the first order in £ we can ignore the difference between £ and 

Given the smallness of £ it is not surprising that chiral 
perturbation theory works remarkably well in describing the 
physical world. An important question in the field is to find 
the range in £ where 1-loop perturbation theory will be valid 
up to a given error say 1% 0, !!• This will help lattice 
QCD calculations to extract reliably the low energy constants. 
Current lattice calculations typically use 2-loop chiral pertur- 
bation theory in the region 0.02 < £ < 0. 1 to fit the data in or- 
der to extract the low energy constants of QCD ll^ [Tol[TlT[T2Tl . 
Is this reasonable? What is the physics that controls the con- 
vergence properties of the chiral expansion? Answers to such 
questions are crucial for future progress. Typically one be- 
lieves that it is the p meson resonance that puts the limit on 
pion masses where chiral perturbation theory will be valid. In 
order to avoid physically important singularities in ix — ir scat- 
tering, it is reasonable to expect M n < M p /2 is necessary for 
chiral perturbation theory to be reliable. Experts believe that 



perhaps one needs at least M n < M p /3 |0]. 

In principle, there is another resonance that can limit the 
convergence of the chiral expansion. This is the so called cr- 
resonance and arises in it — tt scattering in a channel with 
vacuum quantum numbers. Recently, the properties of this 
resonance in the physical world were estimated using experi- 
mental input, dispersion theory and chiral perturbation theory. 
It was estimated that M a ~ 440MeV and ~ 544MeV 
O. This makes it a very broad and perhaps not so inter- 
esting resonance in the context of the convergence of chiral 
perturbation theory. On the other hand, in lattice QCD, as the 
pion masses increase, this resonance could become sharper 
and a could become a stable physical particle, a bound state 
of two pions. However, it could also remain an uninteresting 
resonance up to much higher pion masses. Recent studies find 
that the properties of the cr-resonance do depend strongly on 
the quark mass ITilfTsll . It is interesting to ask if this depen- 
dence can affect the chiral expansion. Although this is a diffi- 
cult question to answer in QCD, it may be possible to explore 
it with simpler models. Here we show that the cr-resonance 
can in principle affect the chiral expansion. In particular we 
demonstrate that a light and weakly interacting a can induce 
an early break down of chiral perturbation theory. 

It is easy to argue that a light cr-resonance can indeed trig- 
ger the breakdown of chiral perturbation theory. Consider a 
non-linear sigma model which contains a coupling T that can 
be tuned such that for T < T c it is in a phase where the global 
symmetry is spontaneously broken and for T > T c it is in a 
symmetric phase. Chiral perturbation theory must be useful in 
describing the low energy properties of the theory in the bro- 
ken phase, but not in the symmetric phase. This means, as T is 
tuned towards T c in the broken phase, chiral perturbation the- 
ory must become poorly convergent. Close to T c , if the phase 
transition is second order, the linear sigma model becomes a 
good description of the physics and in that model, as we will 
see later, the breakdown of chiral perturbation theory can be 
traced to the fact that M^/F^ becomes small. Note that at 
the critical point the sigma and the pions become degenerate 
and chiral symmetry is completely restored. Although the sce- 
nario that a light cr-resonance affects the convergence of chiral 
perturbation theory is perhaps known to the experts, we do not 
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know of any previous work which demonstrates this explic- 
itly. This is the main motivation for our current work. Here 
we study a QCD-like lattice field theory model which has the 
same symmetries as two-flavor QCD. Hence SU(2) x SU(2) 
chiral perturbation theory is applicable. Our model also con- 
tains a parameter equivalent to the coupling T of the non- 
linear sigma model discussed above. We tune this coupling 
to be close to the critical point and hence know that our model 
contains a light sigma resonance although we do not know its 
exact properties a priori. We then find evidence that indeed 
chiral perturbation theory breaks down when M n > M a /3 is 
roughly satisfied. 

Our model involves two flavors of staggered fermions inter- 
acting strongly with abelian gauge fields. We recently devel- 
oped an efficient cluster algorithm for this model and studied 
it in the e-regime 11611 . Here we will focus on the p-regime. 
The action of the model is given by 
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where x denotes a lattice site on a 4 + 1 dimensional hyper- 
cubic lattice L t x L 4 . Here L 4 is the usual Euclidean space- 
time box while L t represents a fictitious temperature direc- 
tion whose role will be discussed below. The two compo- 
nent Grassmann fields, ip x and ij) x , represent the two quark 
(it, d) flavors of mass m, and (p^ <x is the compact £7(1) gauge 
field through which the quarks interact. Here /i = 1, 2, .., 5 
runs over the 4 + 1 directions. The /j, = 1 direction will de- 
note the fictitious temperature direction, while the remaining 
directions represent Euclidean space-time. The usual stag- 
gered fermion phase factors rj^^ x obey the relations: rjf x — T 
and T]f x = l for i = 2,3,4,5. The parameter T controls 
the fictitious temperature. The four fermion coupling c sets 
the strength of the anomaly. As explained in [16], the above 
model has the same symmetries as Nf = 2 QCD, i.e.., when 
to = 0, the action exhibits a global SUl(2)xSUr(2) symme- 
try, which is explicitly broken down to SUv{2) when to ^ 0. 
In this work we fix L t — 2 and c = 0.3. For these parame- 
ters the temperature T can be tuned so that the model is in a 
spontaneously broken phase for T < T c or in the symmetric 
phase for T > T c , where T c = 1.73779(4) was determined 
in the earlier work lfl6ll . Since the phase transition is second 
order, close to T c the pion decay constant in the chiral limit 
F is small in lattice units. This reduces the lattice artifacts 
in our model. Further, tuning T close to T c also makes the 
CT-resonance light as discussed above. For these reasons, we 
chose to fix T = 1.7 in this work. 

We focus on three observables: The vector current suscep- 
tibility Y v , the chiral current susceptibility Y c and the chiral 
condensate susceptibility \a- The current susceptibilities are 
defined as 



Y v c — 



dL d\12 ( 12 j m c 
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where (x) and J° (x) denote one of the components of the 
vector and the chiral current respectively. The condensate sus- 
ceptibility is defined as 
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For a detailed discussion of our algorithm and observables, we 
refer the reader to lfl6ll . 

The behavior of these observables for large L and small m 
is governed by chiral perturbation theory which is described 
by the Euclidean chiral Lagrangian density 



C=^-tv(dM^daU)-^tr(U + U^ 



(4) 



where F is the chiral pion decay constant, S is the chiral con- 
densate and U £ SU(2) is the pion field. Using this La- 
grangian, the finite size scaling formulas for many quantities 
have been found in in the literature 03 EH EH El El- The 
pred ictions for Y c , Y v and \a in the p-regime can be found in 
MM: 
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where M OT is the pion mass, F v is the pion decay constant 
and (qq) is the chiral condensate at a given quark mass to. 
The function <fi arises due to pions constrained to be inside a 
periodic box and is given by 
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where K\ is a Bessel function of the second kind and n — 
n\ + n\ + n\ + n\. 

We have varied the quark mass in the interval 0.0002 < 
to < 0.01 for lattices in the range 12 < L < 32. Our data fits 
well to the above predictions of chiral perturbation theory for 
0.0002 < to < 0.0035. The detailed results are summarized 
in Table U Thus, we are able to extract F„, M n and (qq) as 
functions of the quark mass. As an illustration, we show the 
data at to = 0.0002 and to = 0.001 in Fig. [T] Note that 
the fit is not as reliable at the lowest mass (to = 0.0002) as 
compared to higher masses. It is possible that our lattices are 
not sufficiently large at this tiny quark mass to allow us to fit 
to 1-loop results. 

At to > 0.002 the fits converge only if we exclude almost 
all the curvature in Y c and \<y- I R particular, we are not sen- 
sitive to the <7i(A) function for these two observables and the 
data fit well even to a constant as shown in Table HTl We illus- 
trate this issue by plotting the data and the fits at m = 0.0035 
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TABLE I: Results from fitting Y V ,Y C , and \ 


CT as a function of L to 
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the finite-size one-loop chiral perturbation theory. The x quoted is 
per degree of freedom 
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0.2536(2) 0.9 
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0.2521(5) 1.8 



TABLE II: Results from fitting Y c and Xa to a constant while Y v is 
fit to one-loop chiral perturbation theory. 
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FIG. 1: Finite size scaling of Y V ,Y C and x<r at m = 0.0002 (top) 
and m = 0.001 (bottom). The solid lines are fits of the data to the 
expected form from chiral perturbation theory. 



and 0.0065 in Fig. [2] Comparing the results from the two 
different fits we see that the error bars for (qq) are underesti- 
mated by a factor of two or three at the higher masses. We find 
that Mjr can be calculated very accurately by a one-parameter 
fit of Y v which may be a useful observation for lattice QCD 
calculations. Interestingly, Y v continues to fit well to the one- 
loop formula even at higher masses, but n (in Eq.|6j could be 
restricted to small values (typically less than 3). 

The quark mass dependence of F„, (qq) and have been 
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CM 



0.4354(3) 0.2329(2) 11.9(3) 19.3(5) 39(3) 
0.4351(5) 0.2331(4) 12.3(5) 18.9(9) 37(3) 



1.1 
1.6 



TABLE III: Results from a combined fit of the data in Table U to 
Eqs.|7] The first row uses data in the range 0.0002 < m < 0.001 
while the second row excludes m = 0.0002 from the fit. 



computed up to 1 -loop in iTTvL [Tell : 



F„ = F 



(qq) = S 



Mi = M 



l-e'log£' + 2£'c F 



l-^'logcT + 3cf Cs 



(7a) 
(7b) 
(7c) 



where cp, cs and cm are higher order low energy constants 
and are usually defined in the literature as Ci = log(A i /47rF). 
We have performed a combined fit of all the values of F^^q) 
and M w quoted in Table [Jin the region 0.0002 < m < 0.001 
to the above three relations. The result is tabulated in the first 
row of the Table [III] We note that the values of F and S 
agree nicely with F = 0.2327(1) and S = 0.4346(2) com- 
puted earlier at m = ifTrjll . Further, in the e-regime we find 
cm + 4cs = 80(6), while in the p-regime (from Table [TTTb 
we see that this number is 87(4). Thus, we confirm that the 
p-regime and the e-regime are described by the same low en- 
ergy constants as expected. This is the first important result of 
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FIG. 2: Finite size scaling of Y V ,Y C and \ a at m = 0.0035 (top) 
and m = 0.0065 (bottom). The solid lines are fits of the data to 
the expected finite size scaling form from chiral perturbation theory 
while dashed lines are fits to a constant. 
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FIG. 3: Rescaled and subtracted quantities Rf„, R(qq) an d Rm 2 
defined in Eqs. l[8j. The solid lines are plots of the fits discussed in 
the text. The dashed lines show the linear region for larger values of 
The "knee" is estimated roughly as the point where the two lines 
cross. 



our work. 

In order to isolate the region where one-loop corrections are 
a good description of the data we define the following rescaled 
and subtracted quantities: 



R F = F*/F-l + Z']0gZ', 



(qq)/X -l + 3?log?/2, 



R M = Ml/M 2 - l-g'log(0/2. 



(8a) 



(8b) 



(8c) 



We use chiral values F = 0.2329 and £ = 0.4354 to compute 
the i?'s and By definition, the R's must be linear in £' in 
the region where one-loop results are valid. In Fig.[3]we plot 
the R's as a function of Assuming errors of 1% or less 
can be tolerated, Fig.[3]shows that the linear region of 1-loop 
chiral perturbation theory occurs roughly when £' < 0.002. 
Interestingly, there is also an approximately linear region for 
£' > 0.006 but with a completely different slope. This is 
shown as the dashed line in Fig. [3] This behavior suggests 



that chiral perturbation theory begins to break down. We will 
argue below that the cr-resonance is responsible for this break 
down. Note that £' sa 0.0035 is the rough location of the 
"knee" that separates the low £' and high £' regions. 

The unnaturally large values of cp, cy, and cm are clearly 
responsible for the break down of the chiral expansion at very 
small values of What is the physics behind these large 
values? It has been argued in the context of the 0(4) linear 
sigma model, that the physics in the sigma channel is directly 
related to these terms. In particular, perturbative calculations 
show that iSIH IH : 



(9a) 
(9b) 

(10) 



c s = log(M iJ /47rF) 



7 8^ 

6 3g R 

. 7 8tt 2 

c M = log (Mr 4ttF) -- + 

3 9r 



where 



Ml = Ml 



9R 
16tt 2 



(3tt%/3 - 13) 
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Here M CT is that mass of the a particle and cjr is the corre- 
sponding renormalized coupling, g^ = M^/2F 2 . We believe 
that in our model the above relations must be valid at least as a 
good approximation because we are close to the critical point 
where g^ is expected to be small and the perturbative 0(4) 
linear sigma model is a good description of the low energy 
physics. Indeed, using cs = 12 we find that M a /F ~ 2 
while using cm = 39 we again find that M a /F ~ 2. The fact 
that these two agree with each other is a clear vindication of 
our belief. Assuming M a /F ~ 2 and setting the scale of our 
lattice with F = 90MeV we estimate M a ~ 180MeV in our 
model. At f - 0.0035 we find that M v ~ 60MeV. Hence, 
we conclude that when > M a /3 chiral perturbation the- 
ory begins to break down and the physics is better described 
by the linear sigma model. This is the second important result 
of our work. 

Can we learn something about QCD from our work? Al- 
though there are many important differences between our 
model and QCD, the main difference is that we have tuned 
our model so that it contains a light and most likely narrow 
(T-resonance. In QCD the a is expected to be heavier and 
broader. We think this is the difference why our low energy 
constants turned out to be much larger than QCD, which in 
turn affected the convergence of chiral perturbation theory. 
Clearly, at the minimum we have learned that the properties 
of the cr-resonance do affect the low energy constants of the 
chiral expansion. These properties do change with the quark 
mass while, by definition, the low energy constants are inde- 



pendent of the quark mass. This suggests that chiral pertur- 
bation theory may become reliable only in the region of the 
quark mass where the properties of the cr-resonance do not 
change much. This is the most important lesson of relevance 
to QCD from our work. The quark mass dependence of the 
a and the p was recently studied in [15]. In particular it was 
found that the coupling of a to the pions changes significantly 
with M w , A rough estimate suggests that M n < 250MeV 
is necessary for the properties of the cr-resonance to become 
stable. It would be interesting if this is also the region where 
chiral perturbation theory becomes a reliable tool. 

In summary, here we have studied a model with the same 
symmetries as Nf = 2 QCD from first principles and have 
shown that chiral perturbation theory is a reliable tool only 
for small quark masses. In particular we learned that the cr- 
resonance can be important in determining the convergence 
of the chiral expansion. Studying the quark mass dependence 
of the cr-resonance should be useful and can in principle be 
done using lattice QCD and may shed light on the region of 
quark masses where 1-loop chiral perturbation theory is valid 
in QCD up to a given error. 

We thank G. Colangelo for useful discussions about the cr- 
resonance. We also thank C.Bernard, S. Diirr, C. Haefeli, 
F.-J. Jiang, H. Leutwyler, T. Mehen, K. Orginos, B. Tiburzi 
and U.-J. Wiese for helpful comments. This work was sup- 
ported in part by the Department of Energy grant DE-FG02- 
05ER41368. 
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